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Abstract 



,-J^ I We present the next-to-next-to-leading order (aOg) corrections to the 



first moment of the polarized virtual photon structure function gi{x, Q^, P^) 



Qh! in the kinematical region A^ <C P^ -C Q^, where — Q^(— P^) is the mass 

r-| \ squared of the probe (target) photon and A is the QCD scale parameter. 

In order to evaluate the three-loop-level photon matrix element of the 

flavor singlet axial current, we resort to the Adler-Bardeen theorem for 

r> \ the axial anomaly and we calculate in effect the two-loop diagrams for 

c^ ' the photon matrix element of the gluon operator. The aal corrections 

are found to be about 3% of the sum of the leading order (a) and the 

next-to-leading order (aa^) contributions, when Q^ = 30 ~ lOOGeV^and 

P^ = 3GeV^, and the number of active quark flavors nj is three to five. 
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Figure 1: Deep inelastic scattering on a virtual photon in the polarized e^ e collider 
experiments. 

1 Introduction 

The investigation of the photon structure has been an active field of research both 
theoretically and experimentally in recent years ^1 121 El IH E] ■ Also there has been 
growing interest in the study of the spin structure of the photon. In particular, 
the first moment of the polarized photon structure function gl has attracted much 
attention in connection with its relevance to the QED and QCD axial anomaly 
III IHl El Uni • The polarized photon structure functions can be measured from two- 
photon processes in the polarized e^e~ collider experiments as shown in Fig. 1, 
where — Q^(— P^) is the mass squared of the probe (target) photon. The photon 
structure functions are defined in the lowest order of the QED coupling constant 
a = e^/47r and, in this paper, they are of order a. 

For a real photon (P^ = 0) target, there exists only one spin-dependent structure 
function 5'7(x, Q"^). The QCD analysis oi g1 was performed in the leading order (LO) 
(the order a) [TT], and in the next-to- leading order (NLO) (the order aag) [121113, 
where as = 5'^/47r is the QCD coupling constant. In the case of a virtual photon 
target (P^ 7^ 0) there appear two spin- dependent structure functions, gj^x, Q^, P^) 
and g2{x, Q^, P^). The former has been investigated up to the NLO by the present 
authors in [T31[Tn|, and also in the second paper of ^3]- In R,efs.[T^IT^ the structure 
function gi{x, Q^, P^) was analyzed in the kinematical region 

A'<P'<g', (1.1) 

where A is the QCD scale parameter. The advantage in studying the virtual photon 
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target in the kineinatical region (|1.1|) is that we can calculate structure functions 
by the perturbative method without any experimental data input ,16j, which is 
contrasted with the case of a real photon target where in the NLO there exist 
nonperturbative pieces IT7]. 

In the present paper we focus on the photon structure function g1, and especially 
on the first moment sum rule of gi{x,Q'^,P'^) in the kinematical region (ll.lj) . We 
present a new result on the next-to-next-to-leading order (NNLO) (the order aa"^) 
corrections to this sum rule. 

The polarized structure function gj of the real photon satisfies a remarkable sum 

rule El El 13 GDI 

f\U^,Q')dx = 0. (1.2) 

Jo 

In particular, applying the Drell-Hearn-Gerasimov sum rule jTH| to the case of a vir- 
tual photon target and using the fact that the photon has zero anomalous magnetic 
moment, the authors of Ref. ^U] showed that the sum rule, (jl.2j) . holds to all orders 
in perturbation theory in both QED and QCD. 

When the target photon becomes off-shell, i.e., P^ ^ 0, the first moment of the 
corresponding photon structure function gi{x, Q^, P^) does not vanish any more. In 
fact, for the case A^ <C P^ <C Q^, the first moment has been calculated up to the 
NLO (aa,), as follows PEI; 

3a 



' dxgl{x,Q\P^) = 

vr 



i:4fi-"-'«^' 



=1 



TX 



^^^^2,2 f^siP') a.(Q2 



+ 0(«a2^, (1.3) 



with Po = ll — 2nf/3 being the one-loop QCD P function. Here asiQ"^) is the QCD 
running coupling constant, and Cj is the electromagnetic charge of the active quark 
(i.e., the massless quark) with flavor i in the unit of proton charge and Uf is the 
number of active quark flavors. The first moment depends only on the average charge 
squared (e^) = {J2i^ 6^)/n/ and the average of the fourth power of the charge (e^) = 
(Sj ^ ^t)/^f- Note that the r.h.s. of Eq. (jl.3|) does not involve any experimental data 
input. This is because the nonperturbative pieces, which take part in the real photon 
target case, can be neglected in the kinematical region ()1.1|) and the whole twist-two 



contributions to gj [x , Q"^ , P'^) can be computed by the perturbative method. It is 
also noted that the first term in the square brackets of the r.h.s. of Eq. p.3|) resuhed 
from the QED triangle anomaly while the second term comes from the QCD triangle 
anomaly [HI CHI • 

2 Theoretical framework based on OPE 

For the analysis of the NNLO (aa^) corrections to the first moment of gi{x, Q^, P^), 
we will apply the framework of the operator product expansion (OPE) supplemented 
by the renormalization group method. First recall that for the OPE of two electro- 
magnetic (and thus gauge-invariant) currents, only gauge-invariant operators need 
to be included with their renormalization basis ^Hj- Since gauge- invariant twist- 
two gluon and photon operators with spin one are absent, we consider only quark 
operators, i.e., the flavor singlet J^^ and nonsinglet J^j^s axial currents, as follows: 

J^s = ^7^75 IV' , J^NS = i^ri^iQlh - (e')l)^ , (2.1) 

where 1 is an nj x n/ unit matrix and Qlf^ is the square of the Uf x nj quark-charge 
matrix so that Tt{QI}^ — (e^)l) = 0. Writing the photon matrix elements of the 
quark currents as 

ilipMril^'Mp)) = 6l[-2z6'^^>Je. {l{p)\\Mf^')Mp)) , t = S,NS (2.2) 

where e* and e^ are the polarization vectors of the target photon with momentum p 
and /i is the renormalization point, then the first moment sum rule of gi{x, Q^, P^) 
is expressed as 

£dxgUx,Q',P') = (7(p)ll^55(/i')||7(p))C75(gV/^',M/i'),«) 

+ {lipW,Nsil^')Mp)) C^siQVl-^',9il^'),a) . (2.3) 

Here Cs and Cns are the coefficient functions corresponding to the currents J^g and 
J^jv55 respectively. Putting it more closely, Cs and Cns are the n = 1 coefficient 
functions which appear in the OPE of two electromagnetic currents. Throughout 
this paper we neglect the effect of quark masses. 



We choose the renormahzation point at fi"^ = P^. For — p^ = P^ ^ A^, we can 
calculate perturbatively the photon matrix elements of the axial currents, which are 
expressed in the form as 

(7(p)||J5.(/x' = P')Mp)) = £^- t = S,NS, (2.4) 

with 

^,=Af+^A'"+(^rA'^'+-. (2.5) 

The leading terms Ag and Aj^^ are connected with the Adler-Bell- Jackiw anomaly [2IIj 
and are already known i^. In Sec. 4 we will show, using the nonrenormalization 
theorem for the triangle anomaly 



|{1) _ aW _ /l(2) 



^5 — ^NS — ^NS — ■ (2-6) 



In fact, the result Ag = Aj^g = has been used in Ref.^l] to obtain the NLO 



(aas) corrections to the sum rule shown in Eq. p.3|) . On the other hand, Ag is 
nonvanishing and will be calculated in Sec. 4. 

The Q^ dependence of the coefficient functions Cs and C^s is governed by the 
renormahzation group equations. The solutions to these equations are given by 



Q{Q'/P'rg{P'),a)=exp 



-^i^\^M9') 



Q(1,^(Q2),«), i = S,NS (2.7) 



'siQ^) (3{g') 

where '~fi{g) is the anomalous dimension of the axial current J^j and (3{g) is the QCD 
/3-function. We expand ■yi{g) in powers of g as 

^^(9) = ll''-£^,+ll'\-£-,? + l?\^,f + 0{g') , ^ = S,NS. (2.8) 

Since the flavor nonsinglet axial current J^^s is conserved in the massless limit, it 
undergoes no renormahzation, and we have 

iPs = 41 = 7^ = • ■ ■ = . (2.9) 

Thus the nonsinglet coefficient function Cns is expressed as 

Cns{QVP^ 9iP'), «) = CNsil, giQ'), a) . (2.10) 



On the other hand, the flavor singlet axial current J^g has a non-vanishing 
anomalous dimension jsig) due to the axial anomaly. To be precise, we know 
75. = at one-loop, but at higher loops we have nonzero 75 , 75 and so on. The 
f3 function is expanded as: 



dg_ 

Then using Eqs. 
exp 



Pig) = -Po- 



9" 



-Pi- 



9^ 



-P2 



9' 



^(^') , ns{9') 



167r2 ^^'(16712)2 ^"(167r2)3 
and ()2.11|) . we obtain up to the order of a1 

7i'V«,(g2) «,(p2 



+ ■■ 



(2.11) 



^ ^ 2/5o ^' 47r 



A-K 



4/?o 



Lf J2) _ (1) Ax f^osmv _ (^A^y 



4(#)Y^^-^^V + 0(afl 



dg , 



8V/3o^ V 47r A-K J ' ^'^'' ■ ^^-^^^ 

We need the information on the j3 function only up to the two-loop level, i.e., j3q 
and /3i, in the above expression. But later for numerical analysis, we will use the 
QCD running coupling constant asiQ"^) where the three- loop /32 is also taken care 
of |2g. 

Finally the flavor singlet and nonsinglet quark coefficient functions, 6*5(1, g{Q^), a) 
and C]\[s{l,g{Q'^),(y), are expanded in power of as(Q^) up to the two-loop level as, 

Cs(l.S(Q\a) = (e^>{l + 4"fliM + Bf(^)^...}. (2.13) 



CNs{h3{Q'),a) 
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T— }■ 



(2.14) 



An ""V 47r 

Then putting Eqs. (|2.4ll2.(ijl . flJ^ . (|2. 12112.1 4 jl into Eq.(j2ini), we obtain the expres- 
sion for the first moment sum rule of f77(^; Q^^ P^) ^P to the NNLO {aal) corrections 
as follows: 
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47r 
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'as{Q' 



(3on^ 47r 



47r 



47r 



47r 



(2.15) 



3 Parameters in the MS scheme 



All the quantities necessary to evaluate the NNLO [aal) corrections to the first 
moment of gi{x, Q'^,P'^) have been calculated in the literature and already known, 
except for the photon matrix element of the flavor singlet axial current at three 
loops. If not otherwise mentioned, all the expressions listed in this section are the 
ones calculated in the modified minimal subtraction (MS) scheme 



The flavor singlet axial current J^g has an anomalous dimension starting at the 
two-loop order due to the axial anomaly. The two-loop J2S] and three-loop 
results are: 



7? 



UCpnf, 
284 



36Cj)nj 



(1) 



lC,n} 



(3.1) 
(3.2) 



with Cp = 3 and Ca = 3 in QCD. In fact, 75. was first calculated in the Pauli-Villars 
regularization |2ZI- The result 7^ was obtained in the MS scheme. 

It is well known that a suitable prescription is required for the renormalization 
of the flavor singlet and nonsinglet axial currents in the framework of dimensional 
regularization. For the flavor singlet case, i.e., in order to obtain 7^ in Ref. |26j. the 



75-niatrix was defined as [2H1 



7/^75 = ^e/.par7''7''7^- 



(3.3) 



In addition, after tlie standard ultraviolet renormalization of the axial current J^g in 
the MS scheme, an extra finite renormalization was introduced so that the one-loop 
character of the operator relation of the axial anomaly in QCD 



(9 T"^ 
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^ /lU^ pK 1 



(3.4) 



where G^^, = d^A^ — d^A'^^ + gf'^^^A'^^A'l is the gluonic field-strength tensor, is 
retained intact to all orders in perturbation theory. 

As for the prescription for the flavor nonsinglet case [221, the 75-matrix defined in 
Eq. (|3.3p is used but this violates the axial Ward identity, which is to be restored by 
an additional finite renormalization for J^^g. Then we have the nullified anomalous 
dimension for the nonsinglet axial current (see Eq. ()2.9p ). 



The flavor singlet coefficient function (^^(l, (^(Q^), a) was calculated up to the 
two-loop level jHOl Ell while the nonsinglet coefficient function Cj^siXi9{Q'^)i(^) up 
to the three- loop level [21] ■ The one- loop [32] and two- loop [SHI JHH (HH] results are 



B 



(1) 



B 



(2) 



r(1) 

Cf 



-3Cf 



y C^ - 23C^ + (8C3 + y )n/ 



B 



(2) 

NS 



c, 






23C^ + An 



(3.5) 
(3.6) 

(3.7) 



Both Bf and 5^ 



where Cs is the Riemann zeta- function ((^3 = 1.202056903- ■ 

were calculated in the MS scheme with the prescription for the renormalization of 

the flavor singlet and nonsinglet axial currents explained above. 



The /3 function has been calculated up to the four-loop level in the MS scheme 
IHK] . For the purpose of this paper we need the /3 function up to the three-loop level 

(see Eq. ()2.11| l) and, in addition to [3q, we get 
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(3.8) 
(3.9) 



4 The Adler-Bardeen theorem and the photon 
matrix element of the axial current 

In this section we show how to calculate the photon matrix elements of the axial 
currents (Eqs . (12.4112.51) ) up to the three- loop level {aa^). We make a full use of the 
Adler-Bardeen (AB) theorem [22] for the axial anomaly both in QED and in QCD. 
With its use, the loop-level of the Feynman diagrams to be evaluated decreases by 
one and the calculation becomes much simpler. Also the use of the AB theorem is 
legitimate from the viewpoint of the renormalization scheme (RS) dependence. We 
already know that we should use the same RS to compute higher-order coefficient 
functions and anomalous dimensions of operators in order to get the RS independent 
predictions. In the present case, the same RS should be employed also for the 
calculation of the higher-order photon matrix elements of the axial currents. Recall 
that the definition of the 75-matrix in Eq. ()3.3|) within the MS scheme leads to an 
additional finite renormalization required for the axial currents. In particular for 
the fiavor singlet axial current, the finite renormalization is fixed so that the AB 
theorem for the axial anomaly may be secured. Using this RS or renormalization 
prescription the three-loop anomalous dimension 7^ of Eq. ()3.2j) and the two-loop 
coefficient function Bg of Eq. (j3.6|) have been calculated. Therefore it is a consistent 
procedure to apply the AB theorem for the axial anomaly to the calculation of the 
higher-order photon matrix elements of the axial currents. 

4.1 The AB theorem for QED 

The leading terms Ag and Aj^g in Eq. (j2.5|) are already known [^ El ■ 

aP = -12n,{e') , ASi = -12n,((e^) - (e^)^) . (4.1) 

We will rederive the above results as an example of how to use the AB theorem for 
the calculation of the photon matrix elements of the axial currents. 
The AB theorem for the axial current J5 = ipl'^'^b'^ ^^ QED gives 

dnJ^ = ^e'^'^'^F^uFp. , (4.2) 



where F^^, = d^Ay — d^A^ is the field-strength tensor of photon. The term on the 
r.h.s. is the Adler-Bell-Jackiw anomaly [201 • Note that this is an operator form and 
holds to all orders in perturbation theory. Let us define 

X{Q\T[d,J^ix)A\xi)A-'{x2)] |0)Amputated (4.3) 

where A^ is a photon field and T is the covariant time-ordered product. Actually, 
the expression of the r.h.s. does not depend on x. According to the AB theorem, 
this is equal to 



d xid X2e 



i{ipi+P2)x-piXi-p2X2} 



(0|T 4e>'''P^d^A,{x)d,A^{x)A\x,)A-{x2) |0)Amputated . 



Now differentiating both sides of Eq. (j4.3p with respect to pi„, and setting pi 
—p2 = —p |Sni, we obtain 

R''^^{-p,p) = f d^xid^X2 e'P^'''-''^^ 

x{0\T[j^{x)A\xi)A^{x2)]\0)Arapnt.ted 

where we have used an identity i{x — xiYdi^J^{x) = —iJ^{x), since 



(4.4) 



(4.5) 



d X 



i{x - xiyd^,J!^{x) 



[surface terms] — i j d'^x (9^,k(x — xiY J^{x) 
-I I d'^xJ^(x) . 



(4.6) 



Note that the r.h.s. of Eq. ()4.5|) is nothing but the photon matrix element of J^ at 
zero momentum transfer. Differentiating Eq. (|4.4j) with respect to p^, and setting 
Pi = —p2 = —p, we find that R'^^'^{—p,p) is also expressed as 

-P,P) = -^4/rf^xirf^X2e'P(^i-^^) 
^-''P^{0\T\A,{x)dpA^{x)A\x,)A^{x2)]\0), 



R"^ 



xe" 



/ Amputated 



(4.7) 



The leading terms Ag and Aj^g are obtained from the evaluation of the r.h.s. 
of Eq. ()4.5|) in the leading order e^. Instead we compute the r.h.s. of Eq. ()4.7p . It is 
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+ ... /sssss\ + /^^\ + 




(two-loops) (three-loops) 

Figure 2: Two- and three-loop QCD corrections to the basic triangle diagram. 



a tree-graph calculation and we obtain to order e 

16^2 

Now we take into account the quark charge factors and the color degrees of freedom 
and make the following replacement: 



R'^'^{-p,p) = (--^ 4) X [-2te'''^%] . {U 



e 



2 , Q \ V„ „N2 Q„2„ / „2 



^Y.(^^ey = 3e\{e^) for J55 , (4.9) 



e2 _ 3^(e2_^g2^)^g^^)2^3g2^^(^g4^_^g2^2) fo^ j^^^_ (^^^q) 



Finally considering the convention to define the photon matrix element as is shown 
in Eq. (j2.2|) . we reach the expressions (J4.H) for A^ and Aj^g. 

To obtain the NLO and NNLO terms Aj^g and Ajy^ corresponding to the flavor 
nonsinglet axial current J^ns^ ^^ resort to the nonrenormalization theorem j22] for 
the triangle anomaly. The theorem says that there are no radiative corrections (in 
this case, QCD corrections) to the basic triangle diagram. The contributions of the 
two-loop and three-loop diagrams such as shown in Fig. |21sum up to zero. Moreover, 
the three-loop diagrams such as shown in Fig. El in which the axial current part is 

(2) 

connected to the photon- vertex part by two gluon lines, do not contribute to A}fg. 
Thus we have 

A^nI = A^nI = 0- (4.11) 

4.2 The AB theorem for QCD 

In order to evaluate the NLO and NNLO terms Ag and Ag for the photon matrix 
element of J^g, we use the AB theorem for the singlet axial current in QCD, which 
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+ 




+ 



Figure 3: The three-loop diagrams of the gluon-photon scattering type. 



reads 






jXV PK 



(4.12) 



where G" is the gluonic field-strength tensor. This is an operator form and holds 
to all orders in perturbation theory. We define the following amplitude: 

x{0\T[d,.L^six)A\xi)A^ix2)]\0) Amputated (4.13) 

where A^ is a photon field. Differentiating both sides of Eq. ()4.13p with respect to 
Pia, and setting pi = —p2 = —p, we obtain 



S'^'^i-P^p) 



x(0|Tfe(a;)A\xi)A^(x2)l|0) 



Amputated 



(4.14) 



which is nothing but the photon matrix element of the singlet axial current J^g at 
zero momentum transfer. We consider the amplitude S'^^'^{—p,p) up to the order 
e^g^, namely, up to the three- loop level. 

The AB theorem ()4.12p tells us that (pi + P2)kS'^'^'^{pi,P2) is also expressed as 



iPl+P2).S-''iPuP2) 

x(0|T 



9 nf 



167r2 2 



e'^^'"=G^,(x)G'^,(x)A\xi)A^(x2)l|0)Amputated • (4.15) 
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The operator e^'^^^G^^.G'^^ is expanded in terms of gluon fields as 

e^'^P'^G^G^p^ = Ae^'^P^'d^AldpAl + Agf^^e^^'P^dpAlAlAl^ , (4.16) 

wliere the quartic term of gluon fields does not appear due to the fact that 

Now we find that, up to the order e^g'^, Eg . (14.151) can be rewritten as 

2 

(pi+P2).5'^'"(pi,P2) = -|-^z/rf4xirf^X2e^«^^+P^)"-P^"^-P^"^^ 

,,^^f^-P-I^Q\T[d^Al{x)dpAl{x)A\xM^{^2)\ |0)Amputatcd , (4.17) 

since photon does not couple to gluon field directly. Again differentiating both sides 
of Eq. ()4.17p with respect to pio-, and setting pi = —p2 = —p, we obtain up to the 
order e'^g'^, 

y,,--^P-(^0\T\A:{x)dpA':^{x)A\xi)A^{x2)] |0)Amputatcd • (4.18) 



The NLO term Ag is derived from the calculation of S"^'^{~p^p) in the order 
e^g'^. Instead of evaluating the r.h.s. of Eq. ()4.14p . which is a two-loop calculation, 
we compute the r.h.s. of Eq. ()4.18|) . Then the calculation reduces to the one-loop 
level but gives null result in this order since gluon and photon do not couple directly. 
Thus we obtain 

Af = . (4.19) 

4.3 Calculation of Af^ 

In order to obtain Ag , we need to calculate S'^'^'^{—p,p) in the order e'^g^. Again 
instead of evaluating the r.h.s. of Eq. ()4.14|) . which is a three-loop calculation, we 
compute the r.h.s. of Eq. ()4.18|) . Then the calculation reduces to the two-loop level. 
The relevant two-loop diagrams are shown in Fig. |31 There are also two counter- 
term diagrams corresponding to the diagrams (a) and (b). Because the quarks are 
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T X T X r 

(a) (b) (c) 

Figure 4: Two-loop diagrams contributing to S'^^'^{—p,p) in the order e'^g^. 



Dirac fermions, there are two directions for the charge flow. The result for each 
diagram is invariant under reversing the charge flow, so each diagram is computed 
only for one arrow direction and then multiplied by a factor of two. We have 
performed calculation in the MS scheme with n = 4 — 2e, and used the latest version 
of the program FoRM jHZI to do the necessary algebra. 

For the moment let us assume that each quark has a charge e. The contributions 
of the diagrams (a), (b) and (c) are 



Diagram(a) = Diagram(b) 



Diagram(c) 
where 



ra^Ufi 
U 2 J 


r a ttsi 

-47T 47T- 


rasHfi 

^TT 2 J 


■ a as- 
-47r47r- 



3C,5-(^)"{l + 16}[-2ze->,] 



(4.20) 



oC^ pb 



2e 



jj^ X2e.4 82 



-p^ 



{- + ^-16C3}[-2^6'^">,],(4.21) 



S' = (47r)^e-^^^ , (4.22) 

and a factor SCp comes from Tr [T'^T'^] with T" being the generator of the color 
group SUci'i)- The contributions of the counter-term diagrams corresponding to 
the diagrams (a) and (b) are the same and each gives, 

vr 2 -I l-47r47r-l ^—p' 



,y[-^-12}[-2ze^'^%] 



(4.23) 
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Summing up all contributions, we obtain in the order of e^g^ 

S-'-i-p,p) = ^[^feCFUj {^ - I6C3 + 12 ln^}[-2^e-'^%] . (4.24) 

We renormalize the photon matrix element S'^^'^{—p,p) at yU^ = —p^ = P^. Also we 
take into account the quark charge factor and replace a with anf{e^). Then we find 



53 



AP = l2{e')Cj.n}{'-^-8Cs}. 



(4.25) 



Apart from the charge and SUci'i) group factors, the result (j4.25p is consistent with 
the QED version of the two-loop calculation by Larin |^ for the matrix element of 
^,,upKQo^^Qa^ between two gluon states (see Eq.(22) of Ref.| 



5 The NNLO (aaj) corrections 

Putting the parameters in Eas. (i;ilH;-i.2ll . (i;-i.5H;-i.7ll . (HP?!) . iKT\\ and K7E^ into Ea. dTTHll . 
we finally obtain the NNLO (aa^) corrections to the first moment of gJix, Q^, P^): 



2 td2\ 



dxg1{x, Q , P 



asm 



TT 



2 /^ 



e, 



2 as{Q' 






a,(p2) a^(Q2 



TT 



TT 



1 ((3, 59 2 w^ 2 



TT 



4/3o ^/5o 



+ 



2n 



A^ 
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«,(P2) «^(g2 



«KP^) «Kg 



TT^ 



-n" 



n 



n 



m^f-^Hic^^m^r^ 
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where we are back to the notations J2i^ ef (= nf{e^)) and J2i^ ^t (= '^f{^'^))- The 
third to eighth terms in the parentheses are the NNLO contributions. The sum rule 
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nf Co CiQ Cip C2Q C2QP C2P 

3 0.22222 -0.12346 -0.098765 -0.34685 0.032922 -0.41201 

4 0.41975 -0.12346 -0.29630 -0.027306 0.011852 -1.1533 

5 0.43210 -0.042405 -0.38969 0.50097 -0.11860 -1.4062 

Table 1: The coefficients in Eq. ()5.2j) for the ffist moment of gi{x, Q^, P^] 



is expressed in the form as 



f dxgJ{x,Q\P') 
Jo 



3a , a,(Q2) a,{P 



= < Co + CiQ h Cip 

IT [ TC n 

+C2q[—^) +C2QP—^ :^ + ^2p[—^) I , (5.2) 

where the coefficients Co, Cig, Cip, C2q, C2qp and C2p depend on the number of the 
active quark flavors, nf. We hst in Tabled the numerical values of the coefficients 
Cj for Uf = 3,4,5. We see that the values of Cj increase or decrease with Uf. We 
also see that the values of C2P are rather large in magnitude compared with other 

(2) 

coefficients. This is due to the effect of A^, the three-loop photon matrix element 
of the flavor singlet axial current J^g. Indeed, without the contribution of ^4^ , the 
numerical values of C2p read as -0.11386, -0.32510 and -0.40405 for n^ = 3, 4 and 5, 
respectively. 

To estimate the sizes of the NLO (aa^) and NNLO (aa^) corrections compared 
to the LO (a) term and the ratio of the NNLO to the sum of the LO and NLO 
contributions, we take, for instance, Q^ = 30 and 100 GeV^, and P^ = l and 3 GeV^, 
although these values of Q^ and P^ may not be large enough to apply the formula 
(|5.1|) or (|5.2|) to the cases of n/ = 4 and 5. The corresponding values of a^ are 
obtained from Rei.^. We get a^iQ^ = 30GeV2) = 0.2048, a^iQ^ = lOOGeV^) = 
0.1762, a,(p2 = iGeV^) = 0.4996, and a,(p2 = 3GeV^) = 0.3211. The results are 
given in Table El Since the numerical values of C2P are rather large in magnitude 
compared with those of other coefficients, the NNLO corrections become large when 
P^ is small but still satisfles the condition ()1.1|) . In fact, when P^ = IGeV^ and 
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QVGeV" 


PVGeV" 


LO 


NLO 


NNLO 


NNLO/(LO+NLO) 


^/ = 


= 3 


30 


1 


1 


-0.107 


-0.0520 


-0.0582 






100 


1 


1 


-0.102 


-0.0505 


-0.0562 






30 


3 




-0.0816 


-0.0250 


-0.0272 






100 


3 




-0.0766 


-0.0234 


-0.0254 


% = 


= 4 


30 


1 




-0.131 


-0.0695 


-0.0800 






100 


1 




-0.129 


-0.0694 


-0.0797 






30 


3 




-0.0913 


-0.0288 


-0.0317 






100 


3 




-0.0886 


-0.0287 


-0.0315 


^/ = 


= 5 


30 


1 




-0.150 


-0.0802 


-0.0944 






100 


1 




-0.149 


-0.0811 


-0.0953 






30 


3 




-0.0986 


-0.0309 


-0.0343 






100 


3 




-0.0977 


-0.0319 


-0.0354 



Table 2: The NLO and NNLO corrections relative to LO and the ratio of the 
NNLO to the sum of the LO and NLO contributions for the first moment of 
gi{x, Q^, P^). The values of the QCD running coupling constant we have used are 
a,(Q2 = 30GeV^) = 0.2048, as{Q^ = lOOGeV^) = 0.1762, a,(p2 = iGeV^) = 0.4996, 
and a,(p2 = 3GeV^) =0.3211. 



Q2 = 30 ^ lOOGeV^ the NNLO corrections amount to 6% {uf = 3), 8% (n/ = 4) 
and 9 ~ 10% {uf = 5) of the sum of the LO and NLO contributions. On the other 
hand, when P^ = 3GeV^ Q^ ^ 30 ~ lOOGeV^ and Uf is three to five, the NNLO 
corrections are found to be about 3% of the sum of the LO and NLO contributions. 

6 Summary 

We have investigated the next-to-next-to-leading order (aal) corrections to the first 
moment of the polarized virtual photon structure function gi{x, Q^, P^) in the kine- 
matical region Q^ ^ P^ ^ A^ in QCD. All the necessary information on the coeffi- 
cient functions and anomalous dimensions corresponding to the quark axial currents 
has been already known, except for the three-loop-level photon matrix element (the 
finite term) of the flavor singlet quark axial current J^g. Instead of evaluating the 
relevant three-loop Feynman diagrams, we resort to the Adler-Bardeen theorem for 
the axial anomaly, Eq. ()4.12j) . Then calculation reduces to the one in the two-loop 
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level. We evaluate in effect the two-loop diagrams for the photon matrix element of 
the gluon operator, the r.h.s. of Eq. ()4.18p . 

The aa'^ corrections are found to be about 3% of the sum of the leading (a) 
and the next-to-leading (aag) contributions, when Q^ = 30 ~ lOOGeV^ and P^ = 
3GeV^, and the number of active quark flavors nr is three to five. 
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